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IntybJuct ion 

'  ''  The  present  study  is  investigating  the  basic  neuromotor  coordination 
mechanisms  invoiced  in  a  rapid  elbow  flexion  movement,  and  in  developing 
mathematical  models  to  explain  the  interaction  of  these  basic  neuromotor  co¬ 
ordinations  with  the  biomechanical  parameters  of  movement  speed.  Speed  of 
movement  is  being  assessed  in  a  biomechanics  mode  via  displacement,  velocity, 
accelerat ion ,  point  of  inflection  between  acceleration  and  deceleration,  and 
the  total  time  of  an  elbow  flexion  movement.  Electromyographic  analysis 
techniques  are  used  to  monitor  the  sequential  timing  of  agonist  and  antagonist 
muscle  activity.  The  mathematical  modelling  effort  incorporates  the  bio¬ 
mechanical  parameters  into  an  interface  with  the  neurophysiological  parameters 
involving  the  central  and  peripheral  nervous  systems,  and  then  extends  the 
interface  to  include  viscoelastic  properties  of  the  muscle,  activation  delays, 
and  neuronal  pools. 

The  experimental  approach  includes  consideration  of  the  neuromotor 
coordination  mechanisms  in  both  loaded  and  unloaded  elbow  flexion  movements; 
changes  in  control  mechanisms  due  to  practice  and  learning  effects;  changes 
in  control  mechanisms  due  to  local  muscular  fatigue  induced  by  isometric 
exercise  in  the  agonist  and  in  the  antagonist  muscle  groups;  feasible 
training  and  practice  regimens  involving  artificial  means  of  enhancing 
beneficial  changes  in  control  mechanisms;  and  development  of  suitable 
mathematical  models  to  explain  in  functional  terms  the  ways  in  which  the 
adaptive  mechanisms  can  account  for  changes  in  basic  coordination  skill  and 
the  breakdown  of  skilled  movement  patterns  due  to  local  muscular  fatigue. 

The  planned  series  of  studies  incorporates  research  protocols  from  at  least 
three  usually  distinct  and  isolated  areas  of  research:  neurophysiology, 
biomechanics,  and  computer  science  and  mathematical  modelling. 


Summarv  of  Part  I : 


The  first  formal  investigation  was  completed  December,  1980  a:.  1 
data  are  presently  being  analyzed.  The  first  study  was  designed  to  secure 
criterion  measure  estimates  essential  for  input  into  preliminary  mathematical 
modelling  attempts  for  an  unloaded  fast  elbow  flexion  movement.  As  the 
early  mathematical  models  developed,  it  was  necessary  to  revise  experimental 
methodologies  and  to  collect  additional  criterion  measures  as  dictated  by 
modelling  needs.  The  first  stage  of  mathematical  modelling  has  been  completed 
and  is  described  in  Part  II  of  this  report. 

Beginning  in  January,  1981  the  second  formal  investigation  was  begun. 
Again  guided  by  the  insights  and  requirements  generated  by  the  mathematical 
modelling,  the  second  study  is  assessing  basic  neuromotor  mechanisms  involved 
in  fast  elbow  flexion:  (a)  in  loaded  and  unloaded  conditions,  (b)  under 
conditions  of  rested  and  induced  fatigue  of  the  agonist  or  the  antagonist 
muscle  groups,  (c)  a  unique  condition  in  which  agonist  and  antagonist  muscle 
groups  are,  simultaneously  fatigued,  and  (dl  in  male  and  female  subjects. 


Part  t-A  Completed  Work 

Both  Principal  Investigators  devoted  full  time  during  the  summer  to  the 
grant  by  (a)  identifying,  reviewing,  and  evaluating  recent  research  publications 
relevant  to  theoretical  aspects  of  the  grant  project,  fb)  planning  the  revised 
research  protocols  in  detailed  methodology,  and  (c)  identifying,  reviewing  and 
evaluating  recent  research  publications  relevant  to  data  analysis  and  computer 
programming.  Two  research  assistants  were  employed  during  the  month  of  August 
to  expedite  these  undertakings.  Preliminary  orders  for  necessary  expendable 
supplies  were  sent  out. 

In  September  additional  graduate  research  assistants  (RA'sl  were  recruited. 

In  Dr.  Kilmer's  case  the  RA  continued  on  with  computer  programming  for  the 
mathematic  il  modelling  aspects  of  the  project.  In  Dr.  Kroll's  case  new  RA’s 
were  indoctrinated  into  research  data  collection  protocols.  The  indoctrination 
phase  included  practice  in  apparatus  usage,  maintenance,  and  calibration  technique 
as  well  as  pilot  testing  research  measurement  schedules  to  insure  reliability 
of  data  collection.  The  indoctrination  phase  lasted  six  weeks  resulting  in 
highly  proficient  research  assistants  capable  of  collecting  reliable  data. 

Beginning  the  end  of  October,  data  collection  for  the  first  project  began. 

The  first  study  was  designed  to  secure  essential  criterion  measure  necessary  for 
input  into  preliminary  mathematical  models.  These  criterion  measures  include: 

...maximum  isometric  strength  of  elbow  flexion  and  extension 

...rate  of  tension  development 

...integrated  electromyographic  activity  at  maximal  and  sub-maximal 
tensions  and  joint  angles 

...agonist-antagonist  contraction  patterns  during  unloaded  sneed- 
of-movement  elbow  flexion 

...effects  of  local  muscular  fatigue  upon  ncuromotor  coordination 


Preliminary  mathematical  modelling  attemots  quickly  showed  that 
additional  criterion  measures  would  he  necessary.  As  a  result  of  the 
modelling  requirements,  upper  limb  volumes  and  biceps  and  triceps  reflex 
time  and  rate  of  muscle  tension  development  were  collected.  Limb  volumes 
are  being  estimated  by  anthropometric  and  water  displacement  techniques. 

The  research  assistant  has  been  trained  in  these  techniques  by  Professor 
Frank  Katch  of  the  department  of  Exercise  Science  who  is  a  recognized 
researcher  in  the  body  composition  area.  A  MEDLARS  literature  search 
revealed  only  a  few  scattered  articles  dealing  with  upper  limb  reflex 
times  with  none  of  the  data  appropriate  for  modelling  requirments.  Suitable 
apparatus  was  designed  and  built  for  upper  human  limb  reflex  testing.  Limb 
volumes  and  reflex  times  have  been  collected  on  all  subjects  tested  to  date. 
At  this  time  12  subjects  have  fully  completed  the  first  phase  of  data  col¬ 
lection  and  the  results  of  being  analysed. 

In  this  first  study  of  Phase  I,  each  of  the  12  subjects  was  tested  over 
ten  experimental  sessions  with  each  session  taking  90-120  minutes.  Eight  of 
these  sessions  involved  speed  of  movement  testing  while  the  other  two  session 
involved  upper  limb  reflex  and  volume  assessments.  Each  of  the  first  four 
experimental  days  included  50  trials  of  rapid  elbow  flexion  speed  of  movement 
to  a  designated  target.  Previous  research  has  shown  that  this  amount  of 
practice  insures  that  all  subjects  become  well-practiced  and  exhibit  con¬ 
sistent,  stabilized  performance  both  in  the  speed  of  movement  as  well  as 
in  the  neuromotor  coordination  mechanisms.  Following  the  establishment  of 
we  1 1 -pract iced  performance,  local  muscular  fatigue  was  induced  by  two  dif- 
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with  a  5-second  intertrial  rest  period--designated  the  5:5  eondition--for 
a  total  of  30  serial  trials  constituted  the  more  intense  fatigue  series. 

A  5-second  maximal  isometric  contraction  with  a  10-second  intertrial  rest 
period- -des ignated  the  5:10  condition--for  a  total  of  30  serial  trials 
constituted  the  less  intense  exercise  regimen. 

These  two  exercise  regimens--the  5:5  and  the  5:10  condit ions- -were 
designated  to  produce  different  levels  of  fatigue  in  the  involved  muscle  group. 
The  5:5  and  the  5:10  exercise  regimens  were  administered  on  separate  occasions 
to  the  agonist  (elbow  flexors)  and  to  the  antagonist  (elbow  extensors)  muscle 
groups  using  a  balanced  order  of  presentation  over  subjects  and  across  test 
sessions  to  minimize  contaminating,  test  order  effects. 

One  of  the  purposes  of  inducing  high  and  low  intensity  local  muscular 
fatigue  in  the  agonist  or  antagonist  muscle  groups  was  to  ascertain  the  role 
of  changes  in  the  peripheral  muscle  state  and  in  peripheral  muscle  afferent 
feedback  to  higher  nervous  system  centers.  Specifically,  the  question  of 
whether  or  not  programmed  central  commands  for  a  fast  ballistic  movement  can 
be  altered  by  the  presence  of  different  levels  of  muscular  fatigue. 

Currently  data  are  being  collected  in  Phase  II  of  the  grant  proposal 
(see  pages  10-11  of  original  proposal).  In  this  investigation  the  same  testing 
schedule  outlined  for  the  Phase  I  study  is  being  replicated  for  elbow  flexion 
speed  of  movement.  However,  to  further  elucidate  the  basic  neuromotor 
coordination  mechanisms  a  load  is  imposed  upon  the  limb  movement  task.  The 
imposition  of  a  load,  of  course,  greatly  affects  the  muscle  activation  time  and 
sequential  firing  of  involved  muscle  groups  and  constitutes  a  more  complex 
movement  task  with  distinctly  more  complex  nervous  system  control  system 
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involvement . 
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In  order  to  insure  identical  loading  t'or  subjects  with  different 
limb  lengths,  an  inertial  loading  technique  is  being  used.  The  natural 
moment  of  inertia  can  be  altered  by  placement  of  light  weights  at  long 
distances  and  heavy  weights  at  short  distances  from  the  fulcrum  point. 

The  protocol  being  employed  utilizes  a  fixed  weight  for  all  subjects  but 
varies  the  distance  from  the  fulcrum  point  to  comply  with  the  condition  of 
imposed  moment  of  inertia.  One  load  condition  being  used  is  equal  to  two 
times  the  natural  moment  of  inertial  while  the  second  load  condition  is 
equal  to  five  times  the  natural  moment  of  inertia. 

In  addition  to  the  two  different  exercise  regimens  of  5:5  and  5:10 
administered  to  the  agonist  and  to  the  antagonist  muscle  groups  on  separate 
occasions,  another  set  of  exercise  regimens  are  also  being  administered. 

Based  upon  mathematical  modelling  considerations,  it  is  of  interest  to 
ascertain  the  effect  of  local  muscular  fatigue  in  both  the  agonist  and  antagonist 
muscle  groups  simultaneously  upon  basic  neuromotor  coordination  mechanisms. 

To  accomplish  such  a  goal  of  producing  fatigue  in  both  agonist  and  antagonist 
muscle  groups,  a  f lexion-to-extens ion  contraction  sequence  was  designed. 

By  using  a  5-second  contraction  of  the  elbow  flexors  followed  immediately 
by  a  5-second  contraction  of  the  elbow  extensors  with  no  intertrial  rest  period, 
fatigue  effects  equivalent  to  5-second  contraction  and  5-second  rest  periods 
can  be  produced  in  both  the  agonist  and  antagonist  muscle  groups.  Similarly, 
if  the  effort  sequence  is  a  5-second  contraction  of  the  elbow  flexors  followed 
by  a  5-second  contraction  of  the  elbow  extensors  followed  by  a  5-second 
rest  period,  the  equivalent  fatigue  effects  of  a  5-second  contraction,  lO-second 
rest  period  exercise  regimen  can  be  produced  in  both  the  agonist  and  antagonist 
muscle  groups  s imul taneous ly. 


These  two  new  exercise  regimens,  designated  the  5/5:0  and  the  5/5:5, 
will  allow  assessment  of  fatigue  effects  in  the  agonist  and  the  antagonist 
muscle  groups  upon  basic  neuromotor  coordination  mechanisms  and  provide  a 
stringent  test  of  the  predictive  power  of  the  mathematical  model  being 
developed.  In  actuality  the  two  new  exercise  regimens  will  produce  slightly 
more  fatigue  than  the  5:5  and  the  5:10  regimens  upon  which  they  were  based 
because  some  degree  of  co-contraction  of  agonist  and  antagonist  muscle 
groups  occurs  in  maximal  isometric  contractions.  Such  co-contraction 
effects,  however,  are  desirable  since  they  will  heighten  the  local  muscular 
fatigue  produced. 

The  Phase  II  study  described  above  will  involve  12  male  and  12  female 
subjects  with  each  subject  attending  12  test  sessions  of  90-120  minutes  each. 
Ten  of  these  sessions  involve  speed  of  movement  and  exercise  conditions  while 
the  other  two  sessions  are  for  upper  limb  reflex  and  volume  testing.  One  of 
the  research  assistants  (Zulma  C.  Garcia!  will  be  using  some  of  the  data  for 
her  Ph.D.  dissertation  in  exercise  science  under  the  chairmanship  of  PI 
Troll.  It  is  anticipated  that  data  collection  will  be  completed  by  May,  I9SI 
with  subsequent  reduction  and  analysis  of  data  taking  place  over  the  summer 
months.  A  Master's  thesis  is  currently  being  completed  which  is  a  pilot  study 
effort  dealing  with  the  effects  of  vibration  upon  speed  of  movement  and  neuro¬ 
motor  coordination  mechanisms.  This  graduate  student,  Marilyn  Tevos,  was  not 
a  research  assistant  but  her  work  was  supported  in  part  by  purchase  of  essenti 
expendable  supplies. 


PART  II 
yodel 1 i na 


Sunmary  of  Part  II: 

We  describe  a  two-compartment  model  of  the  neuromuscular  systen 
involved  in  the  voluntary  fast  arm  movement  to  a  target  aiscussec  by  Kr til 
in  Part  I.  The  first  cccartmert  accepts  averaged  biceps  and  triceps  £.VG 
signals  as  inputs,  and  models  the  arm's  musculo-s.keletal  response  by 
producing  elbow  angular  position  .g(t)  and  velocity  ?(t)  over  the  correspcr 
ing  movement  time.  This  is  called  our  E/Q/*(t)  model ,  and  its  defining 
equation  is  of  the  form 

i(t)  =  (extensor  torque  =  Q-)  -  (flexor  tcrque  =  Qp). 

A  detailed  version  of  this  equation  is  being  computer-simulated  with  the 
PASCAL  programme  listed  in  Appendix  A. 

Our  second  compartment,  called  our  cjm/cont/E  model,  accepts  volition 
commands  as  inputs,  and  models  tne  nervous  system's  response  by  producing 
smoothed  flexor  and  extensor  EVG  signals  to  feed  into  the  E/C/;  model.  F; 
1 1-4  gives  a  schematic  representation  of  the  ccm/cont/E  model,  whic.r  is 
explained  in  detail  in  the  report. 

The  purpose  of  our  two  models  is  to  separate  the  various  control 
functions  in  the  ballistic  arm-movement  system  well  enough  to  understand 
where  adaptation  occurs  as  speed,  precision,  and  coordination  of  movement 
improve  with  practice. 

This  is  the  first  publication  that  has  arisen  from  work  or  th;s 
contract . 
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1 1 -A  Completed  Work 

Since  last  June  we  nave  completed  an  initial  specification  of  two 
models  that  will  help  us  to  better  understand  and  guide  the  direction  of 
Kroll's  experiments.  The  first  model  accepts  smoothed  EVG  inputs  to  flexor 
and  extensor  muscles  and  calculates  arm  angle  dynamics  fcr  the  fast  arm  move¬ 
ment  described  above  by  K-oll .  This  model  denoted  the  E/Q/s't)  motel 
(short  for  EMG/elbcw  torcue  =  Q/elbow  angular-  velocity  dQ/dt  =  ?}.  The 
second  model  feeds  flexor  and  extensor  smoothed  EMG  outputs  into  the  first 
model,  and  receives  its  inout  commands  frcm  a  volitional  fast  arm  movement 
center  in  the  brain,  presumably  via  the  pyramidal  tract  frcm  the  cerebellum 
to  the  brainstem  and  spinal  cord  (cf.  Miles  and  Evarts,  1979).  This  model 
is  denoted  the  Com/Cont/E  model  (short  for  command/neural  control/EMG) .  The 
appropriate  connection  of  our  two  models  provides  a  complete  overall  model 
for  the  fast  arm  movement  as  c riven  by  volitional  commands. 

Below  we  describe  the  E/Q/j  model  first  and  the  ccm/cont/E  mocel 
second.  Section  1 1 -8  then  proposes  further  studies  fcr  next  year. 


11  The  E/Q/t  model 


This  model  arose  out  of  a  desire  to  fcr-„iate  a  quantitative  acco^rt 
of  the  E/Q/s  part  of  the  arm  movement  to  target  cescrioed  by  Kro’l.  Af:e 
reviewing  the  relevant  literature,  we  rejected  possiole  model  formulat'd” 
at  the  molecular  biophysical  level  as  too  complex  (cf.  Hatze,  19~S,  1980, 
and  Dijkatra  et  al . ,  1973a,  1973b).  In  addition,  for  our  movement,  trey 
would  be  too  difficult  to  specify  quantitatively.  We  also  rejected  motel 
formulations  at  the  spring-mass  postural  level  (cf.  Saxitt,  1580)  because 
they  would  be  too  steady-state  oriente-a  to  allow  an  analysis  of  c-ur  la-'ce 
transient  arm  motion.  Other  control  engineering  mode's  (cf.,  e.g.,  VcP-ue 
seemed  too  general  or  too  coarse-grained  for  our  purpose. 

Figure  1 1 - 1  gives  a  schematic  representation  of  tne  framework  for  th 
E / Q / a>  model  we  chose.  Its  corresponding  tco-level  analysis  formula  is 

vS  =  0E  -  Qc  ,  (i 

where  each  dct  over  $  indicates  a  time  derivative,  and  Qc  and  QF  are  tne 
respective  torques  about  the  elbow  as  exerted  by  the  extensor  and  flexor 
groups  of  muscles  (all  the  muscles  in  each  group  ope-ate  in  approximate 
temporal  and  tensional  unison  throughout  cur  movement  (bagasse,  1975)). 
is  the  arm's  moment  of  inertia  for  rotations  about  tne  fixed  elbow  jc'nt 
shown  in  Fig.  1 1 - 1  (ever  a  i  change  from  6C°  to  about  -30’). 

At  the  second  level  of  analytical  refinement,  we  let 

Qc  *  •'y c  •  'J'tc  ■*'  tc  (2 

where  ;c  is  the  torque  due  to  viscoelastic  braking  of  the  movement 
(Lestienne,  1979),  and  irTr  is  the  maximum  (over  -90°  Q  9u°(  isometric 

torque  produceable  by  the  extensor,  given  the  smoothed  extensor  E.VS  envel 
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Em6_,  recorded  over  the  last  second  (250  ms,  or  several  seconds  would  serve 
as  well,  but  we  need  a  standard).  To  find  we  accept  the  claim  of 

Messier,  et  al . ,  1971,  that  under  static  loads,  ‘the  averaged  EMG  is 
directly  proportional  to  muscle  tension  and  the  constant  of  proportionality 
is  independent  of  muscle  length."  Admittedly,  the  force  properties  of 
muscles  undergoing  rapid  active  or  passive  stretch  or  compression  are  very 
different  than  when  these  muscles  are  being  actively  or  passively  subjected 
to  static  loads.  But  from  what  little  is  known,  the  best  we  can  do  to 
account  for  this  is  to  include  the  factor  appearing  in  (2).  is  a 

fraction  between  0  and  1  that  inserts  the  force-velocity  relationship 
suggested  by  Thorstensson  et  al . ,  1S76,  and  Perrine  and  Edgerton,  1978,  for 
the  human  knee,  and  by  Lagasse's  1975  observations  on  forearm  accelerations 
during  maximally  fast  movements  (quite  like  our  experimental  ones  only 
untargeted) . 

Our  _  function  is 

F  •  = 

*VE^  1  *  1000  ^ 

L 

where  t  is  in  degrees/sec,  and  the  oouble-cusped  brackets  mean 

x  if  x  >0 
0  otherwise 
For  the  flexor  side  we  let 

qF  =  *VF  •  '•’IF 

where  is  similar  to  and  i>„c  =  Equating  to  for  cur 

movement  is  inaccurate  to  some  degree  because  the  flexor  unce'-goes  active 
contraction  and/or  passive  compression  whereas  tne  extensor  unce’'5ces 
passive  or  active  stretcn  (actually,  the  joint  must  oe  stabilised  sc  tue,"e 
are  no  purely  passive  muscle  states,  only  aoproximate  ones). 


Our  method,  then,  reduces  to  choosing  functions  in  (2)  that  best  fit 
the  averaged  data  before,  during,  and  after  enough  trials  by  each  subject 
(hundreds!)  to  enable  him  to  achieve  maximum  speed.  The  resulting  functions 
should  tell  us  some  valuable  quantitative  things  about  the  muscle  properties 
involved. 

At  the  last  level  of  analytical  refinement,  we  let 


*IF(t)  =  ApCOS?(t ] 


kpEMGp(i)hp(t-r)  dx. 


t-u-5 


(5) 


where  ApCOS$(t)  converts  flexor  force  to  torque  according  to  the  geometry  of 
Fig.  1 1 -1 ;  IMp(x)  is  the  smoothed  EMGp  envelope  at  time  t;  kp  is  a  constant 
of  proportionality  that  we  assume  converts  EMGp(x)  into  flexor  force  at  x 
(Messier,  et  al.,  1971);  and  hp(t-x)  is  the  normalized  flexor  twitch 
response  at  time  t  to  a  unit  kpEMGp  impulse  that  arrived  at  time  x. 

Our  hp(t-x)  is  defined  in  Fig.  1 1 -2 ,  where  we  equate  hp(t-x)  with  its 
extensor  counterpart  (for  want  of  knowledge  to  do  otherwise.  This  could  be 
incorrect  by  perhaps  10%  or  so,  especially  after  many  practice  trials).  The 
shape  of  our  h  function  is  a  ^-truncation  (for  economy  in  computing)  of 
Oijkstra  et  al . ' s  1973a)  twitch  response  function  for  biceps  bracnii  muscle. 
But  whereas  their  =  .015  and  =  .15,  we  began  with  3^  =  .05  and  =  .3. 
As  we  continue  to  fit  the  incoming  data,  our  3 ' s  will  doubtless  change  in 
the  direction  of  Oijkstra  et  al.'s. 

Equation  (5)  assumes  linearity  in  the  accumulation  of  twitch  resconses. 
This  is  probably  somewhat  incorrect  for  our  mcvemert.  Yet  Sakitt,  19SC, 
makes  the  same  linearity  assumes  ion  for  an  arm  mcvemert 


s imil ar  to  ours  and 
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2.  Definition  of  n,(t-r)  =  h.(t-r)  in  equation  (5).  See  text  for 
explanation. 
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obtains  good  results.  In  the  future,  however,  we  plan  to  try  letting  h(t-r) 
vary  as  a  function  of 

t-T 

f  -a(t-T-a) 

1  EMG( a)e  do,  (6) 

t-T-v 


the  exponentially-weighted  EMG  input  from  t-T  back  to  t-t-v  in  time.  This 
seems  like  the  most  expeditious  way  to  introduce  an  appropriate  nonlinearity 
for  the  accumulation  of  twitch  responses  into  equation  (5). 

Letting  Apkp  =  Cp  in  (5)  we  get 


tIF(t)  =  CF  cos*(t) 


t-u 

EMGp(T)hp(t-T)dT 

V* 

t-u-e 


(7) 


Oefining  ^p  similarly  and  neglecting  tp  in  (2)  (as  does  Sakitt,  1981.  cf. 
also  Lestienne,  1979),  our  final  refined  version  of  equation  (1)  becomes 


t-u 


?(t)  =  <ovcos*(t) 


C 


F 


T 


J  £MSp(T)hr(t-x)dt  - 
t-u-5 

t-u  _ 

' 

t)hp( t-T)dt 

J 

t-u-C 


(8) 


with  the  constraining  assumptions  given  above.  To  a  first  approx inat icn  we 
can  obtain  Cp  and  Cp  in  (8)  by  replacing  Ji'  by  the  isometric  torcue  found 


first  in  the  flexor  and  then  in  the  extensor  directions  for  each  subject. 
There  are  various  good  ways  of  measuring  each  subject's  J  (e.g.,  Hatze, 
1980,  refers  to  one  he  has  developed),  but  we  have  not  yet  selected  one. 
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We  now  discuss  how  the  EMti  signals  in  (8)  are  obtained.  Figure  I 1-3 
shows  an  example  test  record.  The  EMGav  signals  are  rectified  and  filtered 
raw  EMG  signals.  EMGav  signals  are  derived  in  an  on-line  computer  which 
"continuously"  leaks  off  a  fraction  of  the  present  integral  of  the  rectified 
EMG  signal.  Denoting  the  latter  EMGr,  EMG  in  Figure  1 1-3,  to  a  good 
approximation,  is  given  by 

t 

EMG(t)  =  J  EMGr(T)e‘a(t‘x)dT  (9) 

0 

Thus  if  EMGr(x)  =  K,  a  constant,  then 

EMG(t)  =  £  [1  -  e‘at]  (10) 

as  t  gets  large.  Also,  if  EMGr ( t)  =  Ct,  C  a  constant,  then 

EMG(t)  =  — t  -  — ?[1  -  e"at]  *  -  t  (11) 

a  c  a 


for  small  t.  Finally,  it  is  easily  shown  that  if  EMGr 
long  time  and  then  EMGf  becomes  0  at  t, 

1  .  .  Ke-«  *  _K 
de 


=  K  as  in  (10)  for  a 


(12) 


for  small  e.  Comparing  (10)  and  (12)  gives  us  a  way  to  determine  a  from 
appropriate  data  samples. 

Equations  (10)  -  (12)  enable  us  to  find  easily  the  feMG(t)  Quantities  in 

A 

(8)  from  EMG(t)  recoras  such  as  shown  in  Fig,  1 1-3.  For  example,  given  the 
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EMGp  signal  in  Fig.  1 1-3 ,  we  form  a  straight-line-segmental  approximation  to 
it  by  connecting  a  to  b,  b  to  c,  etc.  If  in  this  approximation,  a  segment 
is  essentially  constant  at  K,  we  let  l-MG  there  equal  Ka.  If  the  segment 
rises  with  slope  C  >  0,  we  let  EMG  rise  correspondingly  with  slope  Co.  If 
the  segment  falls  from  K  with  slope  -K,  we  let  IRT  there  equal  0.  Thus  our 
£MG(t)  equals  a  times  the  straight-line-segmental  approximation  to  EMG 
except  where  that  approximation  falls  off.  The  fall  off  intervals  are  easy 
to  handle  because  the  EMG  signals  for  our  movement  are  either  on  at  high 
intensity  or  nearly  off  over  most  of  every  trial. 

There  is,  however,  one  adjustment  to  the  "”MG(t)  signal  derived  above 
that  we  need  to  make.  In  order  to  compensate  for  EMG  inaccuracies  during 
EMG  falloff  phases  and  intrinsic  EMG  distortions  of  nervous  excitation  to 
muscle  at  low  signal  levels,  we  shall  augment  our  EMG  signal  with  a  *(t) 
postural  feedback  term.  This  is  justified  by  assumption  4  in  the  next 
section. 

To  make  the  adjustment,  we  denote  the  previous  cMG  by  tMG* ,  ana  add  a 
feedback  term  denoted  IMG**.  This  gives  us 

lM£(t)  =  lMG*(t)  +  lMG**(t;  (13) 

on  both  the  flexor  and  extensor  sides.  We  define: 

fH5J*(t)  =  PF1  j\(t-20)  +  PF2i(t-20)  +  PF3i*(t-20)j  G ( t )  (14) 

where  the  double-cusped  brackets  mean  the  same  as  in  (4);  where  the  quantity 
inside  these  brackets  is  *(t)  as  predicted  from  *,  «,  and  V  20  ms  earlier 
(20  ms  being  the  total  delay  around  the  loop  from  y  receptor  to  a  stimulation 
of  muscle);  where  G(t)  =  0  until  t  first  reaches  1  degrees,  and  G(t)  =  1  there- 


after;  and  where 

the  Pp^  are  constants  (now  set  at  ?F,  = 

rF2  F2 

—  recal 1  that  ? 

is  in  degrees/ms  and  V  is  in  degrees/ms 

).  Similarly, 

c  /  •  v 

a 

TMGT(t)  = 

P£1  -U(t-20)  +  P£?*(t-20)  ♦  P E 3 V i'  t-20) 

)  G(t)  (15) 

*  L  \  “  / 

J 

L 
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where  =  Pp. . 

Final ly,  to  keep  (14)  and  (15)  from  causing  everlasting  $  oscillations, 
we  change  (1)  to 

05(t)  =  q£(t)  -  qF(t)  -  c;i(t)  (16) 


I 


where  the  C  viscoelastic  damping  factor  (cf.  lestienne,  1979)  is  set  “near" 
the  level  that  just  prevents  a  second  $  overshoot  of  the  0  degree  target 
position. 

We  now  consider  adaptation. 

As  our  subjects  accumulate  more  experience,  we  know  from  past  results 
(Lagasse,  1975,  e.g.)  that  their  movement  times  will  decrease.  If  the 
adaptations  enabling  this  occur  partly  in  the  muscles,  this  should  require 
us  to  modify  the  h  functions  in  (8).  Since  in  previous  pilot  experiments 
our  flexor  movement  adapted  several  times  faster  and  further  than  the 
opposite  extensor  movement,  it  could  happen  that  only  hp  in  (8)  adapts  with 
practice.  It  also  might  develop  that  all  adaptation  occurs  in  the  nervous 
system  (cf.  the  com/cont/E  model,  next  Section). 

The  crucial  feature  of  our  E/Q/t  model  is  that  eventually  it  should 
locate  and  characterize  the  adaptive  nonlinear  effects  involved  when  speed 
of  movement  is  increased  with  practice.  This  should  give  us  a  mechanical 
perspective  on  how  to  train  for  speed,  precision,  orace  and  sklll--in 
short,  for  well  coordinated  fitness. 

We  have  included  as  Appendix  A  a  listing  of  the  PASCAL  programme 
written  by  Valerie  Congdon  to  compute  *(t)  and  *(t),  as  outlined  above, 
using  a  file  containing  the  EMGp( t)  and  EMGp ( t )  data  for  one  movement.  The 
program  is  well  documented,  hierarchically  organized  into  procedural  and 
functional  levels  identical  to  the  analysis  levels  in  this  report,  and 
specified  so  as  to  be  easy  to  change.  Ms.  Congdon,  who  will  leave  I'Mass 
*'n  January,  1981,  has  done  an  exemplary  job. 


By  the  end  of  February  1981  the  modelling  project's  RA  will  change  from 
Ms.  V.  Congdon,  who  is  leaving  school  to  accept  a  job,  to  R.  Pelosi,  a  Ph.O. 
student  in  the  IE/Operations  Research  Department.  Ms.  Congdon  will  complete 
the  EMG-torque-movement  simulation  programme  before  leaving.  Her  work  has 
been  of  the  highest  professional  caliber  (cf.  Appendix  A,  a  listing  of  our 
PASCAL  simulation  programme).  Mr.  Pelosi  will  need  much  of  the  spring  term 
after  March  1  to  train  for  the  continuation  of  Congdon's  work. 


2. _ The  Com/Ccnt/E  Model 


This  model  concerns  the  functional  organization  of  the  neural  control 
system  for  the  £/Q/< p  model  described  above.  It  is  an  attempt  to  make 
precise  only  those  conclusions  that  have  been  formulated  by  several  experi¬ 
menters  as  necessary  consequences  of  their  tests  on  fast  arm  movements 
involving  the  elbow  and  sometimes  one  other  joint.  The  com/cont/E  model  is 
amenable  to  few  significant  neuroanatomical  interpretations.  Eventually  we 
hope  to  obtain  a  good  structural  model  that  does  not  suffer  from  this 
defect,  but  that  may  be  a  couple  of  years  away. 

Our  specification  of  the  com/cont/E  model  arose  out  of  the  following 
Assumptions 

1.  For  about  the  first  100  ms  at  least,  our  fast  arm  movement  is  entirely 
feedforward  (is  independent  of  any  feedback),  and  thus  is  unaffected  by 
stretch  reflex  or  tendon  afferents  (cf.  Wadman,  et  al . ,  1979;  Desmedt 
and  Godaux,  1978). 

2.  Selecting  the  levels  Ip  and  of  flexor  and  extensor  EMG  activity 
require  the  first  volitional  decisions.  (As  noted  by  Kroll  above,  all 
of  the  flexor  muscles  act  approx imately  in  temporal  and  tensional 
unison  with  the  biceps  during  our  arm  movement.  Likewise  for  the 
triceps  and  the  extensor  muscles.)  These  levels  are  selected  to  set 
the  speed  and  amplitude,  M,  of  the  movement  with  allowances  for  expected 
inertial  loading  of  the  limb  and  estimated  fatigue  of  the  muscles  being 
activated  (cf.  Lestienne,  et  al . ,  1979;  lestienne,  1979;  Wadman,  et 

al . ,  1979;  Wadman,  et  al . ,  1980). 

3.  The  next  movement  decision  is  to  set  the  relative  timing  of  the  first 
flexor  and  extensor  EYG  bursts.  This  determines  the  movement's  comtina- 
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t ion  of  magnitude  and  speed,  again  allowing  for  the  expected  loading 
and  muscle  fatigue  (cf.  references  in  2  above).  Force-velocity  and 
force-length  effects  are  probably  also  included  in  this  reckoning 
(Thorstensson,  et  al.,  1976;  Perrine  and  Edgerton,  1978),  as  well  as 
viscoelastic  braking  (Lestienne,  1979;  Maughan  and  Godt,  1979). 

4.  The  flexor  and  extensor  muscle  tension  levels  at  the  end  of  the  movement 
epoch  when  the  postural  maintenance  phase  is  entered  are  variable,  but 
the  ratio  R  between  these  two  levels,  (as  reflected  by  their  leaky- 
integrated  EMG  levels)  is  a  fixed  postural  command,  and  is  not  affected 
by  the  starting  point,  direction,  amplitude,  or  speed  of  movement 
(Lestienne,  et  al . ,  1979;  Sakitt,  1980).  Presumably  the  EMG  levels 
themselves  at  the  end  of  the  movement  epoch  are  quickly  adjusted  if  need 
be  (how  fast?)  to  obtain  the  desired  postural  stretch  reflex  stiffness 
(Houk,  1976).  Since  the  biceps  tendon  reflex  loop  delay  is  about  35 
ms,  the  stretch  reflex  loop  delay  is  about  50  ms  (Marsden,  et  al . ,  1976) 
and  our  actual  arm  movement  duration  is  about  300  ms,  the  presently 

.  unknown  speed  with  which  the  set  point  of  the  stretch  reflex  is  reset 
and  the  reflex  subsequently  engaged  will  be  imortant  to  discover. 

5.  In  males  at  least,  continued  practice  of  a  fast  arm  movement  reduces 
execution  time  and  alters  the  movement's  EMG  signature  (Lagasse,  1975). 

In  the  light  of  the  above  assumptions,  we  now  give  a  point-by-point 

description  of  the  com/cont/E  model  (Figure  II-4): 

1.  At  the  top  of  Figure  II-4,  R,  l,  M,  ILF,  and  FF  are  commands  issued 
from  the  volitional  decision  center  for  our  fast  arm  movement  (this 
"center"  might  be  distributed  over  a  large  part  of  the  brain).  ,R  is  the 
ratio  described  in  assumption  4  above;  L  times  gc,  L.g c,  and  L.or  are  tn 
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levels  Lc  and  L,-  described  in  2  above;  M  is  defined  in  2;  ILF,  the 
inertial  loading  factor,  and  FF,  the  fatigue  factor,  referred  to  in  2 
and  3  are  appropr iately  scaled  for  multiplication  into  the  subject's 
moment  of  inertia,  J,  for  the  arm  movement  in  question.  Thus  ILF  >_  1, 
with  IFL  *  1  when  no  external  load  is  added;  and  0  <  FF  _<  1,  with  FF  =  0 
interpreted  as  complete  paralysis  and  FF  =  1  as  zero  fatigue  (defined  in 
terms  of  force  development  capability). 

2.  Np  and  Np  are  loosely  interpretable  respectively  as  flexor  and  extensor 
motorneuronal  pools  that  drive  either  the  flexor  and  extensor  arm 
musculature  for  our  movement  or  just  the  biceps  and  triceps. 

3.  At  t  =  0  the  flexor  EMG  is  activated,  t*  is  the  time  at  the  end  of  the 
first  extensor  EMG  burst.  The  time  courses  of  R  and  L  with  respect  to 
t*  are  hypothesized  as  shown  at  the  top  of  the  figure. 

4.  At  t  =  0,  L  turns  NR  on,  which  then  turns  Np  on  (see  6  below),  all  with 
delay  less  than  de.  N^'s  output  firing  rate  rises  and  falls  somewhat 
sporad ical ly,  but  peaks  at  the  L.ge  level.  The  stretch  reflex  pathway, 
SPp,  into  Np  is  intially  inactive  as  noted  in  assumption  1  above,  bp's 
output  equals  the  algebraic  sum  of  its  inputs.  Np  keeps  Np  shut  off 
over  the  reciprocal  inhibitory  pathway. 

5.  At  time  0,  a  shot  out  of  Ng  turns  Np  off  to  end  the  first  flexor  EMG 
burst,  c  and  k  in  the  formula  for  0  are  scaling  constants.  Since  L 
turns  Nn  on,  the  value  of  c  would  be  0  if  the  flexor  muscle  M_  were 

U  r 

perfectly  linear  in  L  and  if  no  co-contraction  of  MR  and  Mr  were 

necessary  to  stabilize  the  elbow  joint  during  movement.  Since  neither 

of  these  conditions  holds,  c  >  0.  After  Np  shuts  off,  N_  is  no  longer 

inhibited  ana  turns  cn.  Then  N~  keeps  Np  shut  off  by  reciprocal 


inhibition. 


6.  Backtracking  a  bit,  only  one  of  NR's  outputs  is  on  at  a  time.  Initially,  l 

I 

L  turns  NR*s  zR  output  on,  which  is  necsssary  to  keep  Np  on.  Later 
N^’s  output  shot  switches  NR's  output  over  to  z^,  which  is  necessary  to 
keep  Nr  on. 

7.  Nr's  input  Tp  provides  NR  with  an  estimate  of  Mp's  tension,  as  computed 
(or  simulated)  by  Nhp.  The  function  hp  is  Nhp’s  estimate  of  M,-'s  twitch 
response  function  hp  (see  Part  II-A-1)  on  the  E/Q/t  model).  In  contrast 
with  the  nervous  system,  we  would  model  hh_'s  output  at  time  t  with 

t 

Np(o)hp(t-o)da.  3y  assumption  4  above,  the  flexor  tendon  organ 
J 
0 

response,  TORp,  is  effectively  unavailable  to  hhp  for  at  least  30-S5  ms. 

Mutatis  mutandis,  NR  1  s  extensor  circuits  are  the  same  as  its  flexor  ones. 

8.  If  our  movement  is  entirely  feedforward,  seme  such  centers  as  Nhe  and 
Nhp  are  necessary  to  control  N_  and  Nr’s  development  of  an  appropr iately 
stiff  arm  posture  at  the  end  of  the  movement  (assumption  4  above).  The 
possibility  that  there  are  Nh_  and  Nhr  circuits  in  the  cereoellum  seems 
to  fit  Eccles1  concept  of  the  cerebellum  as  a  control  reference  in  fast 
voluntary  movements  (Eccles,  1577).  Lagasse,  1975 ' s  results  showing 
that  EMG  temporal  organizations  for  fast  human  limb  movements  don't 
change  with  fatigue  support  the  idea  that  such  movements  are  entirely 

| 

feedforward.  , 

i 

9.  'when  NR‘s  Tc  input  reaches  the  desired  final  level,  NR’s  output  switches 

! ; 

from  zr  to  ip  in  case  Tp  is  not  up  to  the  desired  final  level,  and  to  t 

a 

z.  in  case  T  is  uo  to  its  final  level.  An  active  i  signals  the  enc  of  ); 

trie  movement  arc  subsequent  '■einstatement  of  the  normal  postural  ma’nte-  li 
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nance  mechanisms.  The  desired  final  flexor  and  extensor  tensions  are 
denoted  Tp  and  7^,  and  are  assumed  to  satisfy  7^R  =  7p  =  L  (assumption  4 
above).  Thus  (Tp/Te)  =  R,  and  the  component  in  Fig.  1 1-4  has  suffi¬ 
cient  inputs  to  compute  all  of  its  outputs. 

10.  If  we  assume  some  background  SPp,  SPp,  and  other  noise,  the  sequence  of 
events  in  9  will  not  always  be  as  simple  as  described.  The  first 

Np  or  Np  EMG  burst  might  sometimes  be  cut  off  prematurely,  forcing 
another  one  or  even  two  small  bursts  later  on  (see  the  L  time  course, 
top  right  in  Figure  I»-4).  The  famous  third  EMG  burst  (e.g.,  'wadman, 
et  al . ,  1980)  night  be  accounted  for  in  this  way  —  an  attractive 
hypothesis  because  Kroll's  experiments  to  date  with  highly  practiced 
subjects  usually  don't  show  a  third  EMG  burst.  The  ability  of  the 
control  components  in  Fig  1 1-4  to  learn  (with  lots  of  practice)  the 
amount  of  viscoelastic  braking  that  actually  occurs,  the  force-velocity 
and  force-length  relationships  of  the  muscles,  and  various  nonlineari¬ 
ties  of  the  system,  offer  other  adaptation  possibilities. 

11.  The#-C>  symbol  in  Fig.  1 1  -4  points  to  likely  adaptation  sites  as  the 
same  movement  is  practiced  over  successive  trials  and  days. 

12.  Ng  and  the  (N^,  Nhp,Nhc)  complex  could  be  combinea  instead  of  function¬ 
ally  partitioned  as  shown,  but  that  would  defeat  the  purpose  of  the 
com/cont/E  model. 
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1 1-8  Proposed  Modelling  Studies  Next  Year 
We  wish  to  proceed  cn  three  fronts: 

1.  The  E/Q/$  model  will  be  adjusted,  and  perhaps  revised  or  augmented,  to 
fit  Kroll's  experimental  results.  His  tests  with  different  loads  and 
fatigue  levels  should  be  very  helpful  to  us  during  this  phase. 

2.  The  Com/Cont/E  model  will  be  scrutinized.  In  particular,  the  NQ  and  Nh 
components  will  be  checked  for  plausibility  against  Kroll's  fatigue 

and  load  variational  test  data.  We  now  believe  that  "efferent-copy" 

Mp.  simulators  such  as  the  functions  actually  exist  (perhaps  at  the 
highly  integrative  thalamic  level  atop  the  brain  stem  reticular  forma¬ 
tion,  or  in  the  cerebellum),  but  we  shall  seek  more  evidence  to  support 
this  idea.  We  also  plan  pilot  tests  of  various  kinds  to  further 
elucidate  the  role  of  stretch  reflexes  at  different  times  during  and 
immediately  following  fast  arm  movements.  To  date ,  Uadman  et  al. 

(1979)  have  shown  that  when  a  subject's  planned  fast  arm  movement  is 
unexpectedly  blocked,  the  pattern  of  EMG  activity  over  at  least  the 
first  100  ms  is  the  same  as  before  blocking.  But  we  know  of  nothing 
else  in  the  literature  on  the  dynamics  of  rapidly  changing  the  set 
points  of  stretch  reflexes. 

3.  We  shall  develop  computer  programs  in  PASCAL  for  solving  the  nonlinear 
delay-differential  equations  that  will  arise  when  the  E/Q/$  model  is 
augmented  to  include  feedback  stretch  reflex  information.  Banks 
(1979,  1980)  has  recently  developed  the  requisite  mathematical  methods, 
but  as  yet  PASCAL  programs  to  implement  them  do  not  exist. 

In  1  to  3  above,  we  shall  be  especially  watchful  of  adaptation  parame¬ 
ters  and  characteristics.  In  year  three  of  the  contract,  we  shall  apply 
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the  results  of  the  first  two  years  to  a  search  for  better  ways  to  measure 
and  train  optimally  coordinated  fast  human  limb  movements,  with  an  emphasis 
on  kinesthetic  learning. 


V  Mote  on  RA  Budget  for  Kilmer 

UMass  now  has  no  limits  on  what  RAs  can  make.  The  going  rate  for  a 
full  RA  to  a  grad  student  in  the  Engineering  School  is  at  least  6000  per 
year,  and  for  Ph.D.  students  this  will  probably  soar  to  over  10,000  per 
year  within  the  next  year  (due  to  the  new  UMass-Industry  Program  headed  by  the 
Dean  of  Engineering).  We  will  try  to  hold  the  line,  however,  at  the 
originally  proposed  5700  per  year. 


Personnel  Change  (for  Kilmer) 

Sy  the  end  of  February  1981  the  modelling  project’s  RA  will  change  from 
Ms.  V.  Congdon,  who  is  leaving  school  to  accept  a  job,  to  R.  Pelosi,  a  Ph.D. 
student  in  the  IE/Operations  Research  Department.  Ms.  Congdon  will  complete 
the  EMG-torque-movement  simulation  programme  before  leaving.  Her  work  has 
been  of  the  highest  professional  caliber  (cf.  Appendix  A,  a  listing  of  our 
PASCAL  simulation  programme).  Mr.  Pelosi  will  need  much  of  the  spring  term 
after  March  1  to  train  for  the  continuation  of  Congdon 's  work. 
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APPENDIX  A 


Contents : 

1)  An  example  test  record  of  a  subject  whose  responses  were 
modelled. 

2)  Model  output  from  the  computer  when  the  example  test  records 
EMG  tracesin  1)  were  used  to  derive  the  input  to  the  model. 

3)  Computer  plot  of  $  in  2)  versus  time.  Drawn  in  for  comparison 
is  *  of  1 ) . 

4)  Documentation  for  and  listina  of  the  PASCAL  program  that 
produced  the  output  of  2) . 


NOTE  BENE: 

The  text  of  the  report  was  prepared  before  the  computer  run  for  this 
appendix  was  made.  Some  late  chances  in  the  PASCAL  program  are  NOT  incor¬ 
porated  into  the  report.  They  include  the  new  values  of  the  model 
Darameters  given  in  the  header  of  2),  where  P-j ,  P^,  and  P3  are  the  P 
parameters  of  equations  (14)  and  (15).  Also,  Cr  was  changed  to  equal 
C  ^  until  $'s  first  return  to  0  degrees  after  its  initial  overshoot  of 
the  target  position,  and  to  equal  C,,  after  that.  This  causes  i  to  be 
heavily  damped  upon  its  first  return  to  target  as  required  by  the  test 
data.  Our  comouter  studies  showed  that  with  the  long  delays  in  the 
stretch  reflex  loop,  nonlinearities  in  the  nerve- to-muscle  pathways 
present  or  not,  the  required  damping  could  only  come  from  a  large  increase 
in  vi coelastic  damping.  Thus  we  are  postulating  a  new  mechanism  on  the 
basis  of  our  model  simulation  results. 
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rcczi  :  £  /  a  /  ph::ct  i  s?-s  /  tciwiie  /  angular  oelccity  ) 


PC9ITICN  ("HD 


v  :  oelcci:.'  *  ice  cphidct  »  icc> 

.  :  INDICATES  30  DEGREE  LINE 

00000  :  INDICATES  0  DECREE  L.NE 


13.000 
20.  COO 
DC.  COO 
40.000 
30.000 
SO. 000 
70.000 
30.000 
30.000 
100.00 
110.00 
120 . 00 
130.00 
140.00 
150.00 
130.00 
170.00 
130.00 
ISO. CO 
200.00 
210.00 
220.00 
220. CO 
240.00 
250. CO 
230.00 
270.00 
230.00 
290.00 
300.00 
310.00 
320.00 
230.00 
340.00 
350.00 
2SC.C0 
270.00 

sec. oo 
sec.  co 

400.00 
410.00 
420.00 
420. OC 
440.00 
450.00 
430.00 
470.00 
430.00 
430.00 
500.00 


00 

CO 

00 

00 

00 

00 

00 

0 

0 

0 

0  0 
0  0 


subject's  recorded! 


NOTE: 


Allomna  Far  start-time  uncertainty  because  oF  unknown  initial  bac«3round 
activity  m  tne  subject.  ana  notin3  that  the  relative  timira  ana  snares  aF  tne 
two  PHI  curves  above  can  be  matched  almost  exactly  by  intracuama  ronlmeantis 
inta  t.ie  n<t)  Functions  in  the  intearals  cf  the  PHIDDCT  equation,  t.ie  t-o  P-iI 
curves  ac eve  compare  very  well.  Note  especially  t^c  tfiinael  .s:*P2*?2s0i  so  t.ie 
stretch  reflex  is  net  enaaaedl  alsc»  tne  effect  oF  tr.e  heavy  viscoelastic 
Ca.npina  First  arrears  m  csth  PHI  curves  upon  t.ieir  return  to  cere  oearses  afte 
t.ieir  initial  overshoots.  In  the  suo-ect's  trace. this  da.#?:.ua  First  appears  a; 
aFterPHIDCT  Fi~st  turns  positive.  Cur  model  test  Fit  tne  suspect's  trace  «ren 
lOC.rs.  so  the  ocvious  suaaestion  is  that  as  seen  as  the  spiral  core  receives  a 
?HIDC7*C  Feecoach.  it  aererates  a  nervous  sianai  to  areatl/  increase  visccelast 
Ca.DPina  in  Me  lime. 


mt  f/S. 


Appendix  A  :  PASCAL  Program  for  E/3/PHI  Model 
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Listing  of  PASCAL  PROGRAM  develpoed  far  EMG-Torque— Movement 
simulation  model.  Written  and  tested  by  Valerie  Congdon  in  col¬ 
laboration  with  Prof.  W.  Kilmer,  Fall  and  Winter  1930,  1981. 


Section  Titles  : 

Notation 

Titles  and  Global  Variables 
Procedures  and  Functions 
The  Main  Program 
Input  and  Output 
Exemplary  Output 


... _  .  _  S 


SYMEOL  IN  TEXT  SYMBOL (S)  IN  PROGRAM 


$  >  9  9  5 

Qe.  Qf 
c.  c. 

s 


EMG| 

EMG* 

EMGe,  EMGf 

EMG  =  EMG**  +  EMG* 

G(t) 

hE(t-r) 

hF(t-T) 

q 

qc 


PHI,PHIDCT,PHIDDGT  ;  Position , Velocity .Acceleration 
GUEE  ,  QUEF  ;  Torque  produced  by  triceps, biceps 
ZETA.CZETA  ;  viscoelastic  damping  term  ,  constant 
PSIV(  dot  )  ;  Force-Velocity  Relation 
PSIK  ...  )  ;  Integral  Factor  of  Torque 
ET  ,  GETEMGC  EMGTRICEPS):  Smoothed  EMG  for  triceps 
EB  .  GETEMG(  EMGEICEPS  ):  Smoothed  EMG  for  biceps 
EMGE, EMGF  ;  Corrected  EMG  variables JEXTENSCR , FLEXOR 
EMG  ;  Function  which  returns  Corrected  EMG  values. 
GOFTAU  ;  equals  O.C  if  PHI  >=  C,  1.C  otherwise 
AITCH (T-TAU , EXTENSOR ) ;  Response  of  EXTENSOR 
AITCH (T-TAU, FLEXOR )  ;  Response  of  FLEXOR 
GUE(  MUSCLE  )  ;  q  computed  for  AITCH  function 
QE  ;  q  in  AITCH  Function  for  EXTENSOR 


QF  ;  q  in  AITCH  function  for 

u. 

UE’  UF 

MU  ,  MUE  ,  MUF 

S  1 

5E*  ’F 

XI  ,  XIE  ,  XIF 

A, 

AE’  af 

A  ,  AE  ,  AF 

k. 

Ll. 

UJ 

X  ,  XE  ,  KF 

-  9 

a,  t* 

TAU  ,  SIGMA  ,  T3TAR 

t 

T  ,  ATTIME  ;  time  variables 

sr 

S1E’  S1F 

BETA  1  ,  BETA  IE  ,  BETA  IF 

Co  * 
<_ 

S2E,  £2f 

BETA 2  ,  EETA2E  ,  BETA2F 

e  ,  cos  x 


EXF(X)  ,  CCS(XRADIANS) 


Notation 


OTHER  SYMBOLS  IN  PROGRAM 


TRAPEZOID  ;  Rule  for  integration 
HTRAP  ;  step  size  in  trapezoid  rule 
LASTCNE  ,  LASTFUN  ,  FUN  ;  Hold  the  last 
integration  value  ,  the  last  value  of  t 
the  function  ,  and  the  current  function  value 
HRUNGE  ;  step  size  in  Runge-Kutta  Rule(PSII) 

INITALL  ;  Initializing  Procedure 

SAVEIT .ONTO  =  ARRAY [ 1 . .70 , 0. .2]  of  REAL 

SAVEPHIS  ;  ARRAY  to  hold  values  of  PHI , PHIDCT , PHIDDCT 

PUSH  ;  Procedure  to  save  PHI  values  in  ARRAY  SAVEFHIS 

READEMGFILE  ;  Read  Input  Procedure 

EMGFILE  ;  Input  file 

CUTFILE  .  TAPEC  ,  TAPE  1  ;  Output  Files 

MM  ,  DD  ,  YY  ;  Mounth  ,  Day  and  Year  for  Subject  ID. 

SUB  ,  AP  ,  SES  ;  Subject  No., AM  or  PM  ,  session 

BL  ,  TR  ;  Block  and  Trial  numbers 

WHICHCNE  =  (EXTENSOR  ,  FLEXOR) 

MUSCLE  is  of  type  WHICHCNE 

REAL2D  =  ARRAYCO. .25, 1 . .2]  of  REAL 

EMGA  ,  EMGTRICEP  ,  EMGEICE?  are  of  type  REAL2D 

VALU  ,  TIME  ;  constants  for  array  posticn 

MAX  ,  MAXENGTIME  ;  hold  value  of  maximum  time  to  run 


TITLES  AND  GLOBAL  VARIABLES 
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<*  THE  FASTMCV  PROGRAM  NAS  WRITTEN,  DESIGNED.  AND  IMPLEMENTED  3Y  *) 
<»  VALERIE  CCNGDCN  ,  M.3.  COMPUTER  AND  INFORMATION  SCIENCE,  >) 
(♦  FCR  THE  RESEARCH  BEING  CCNDL'CTED  BY  DPS.  WILLIAM  KILMER  Or  THE  >) 
<*  ELECTRICAL  AND  COMPUTER  ENGINEERING  DEPARTMENT  AND  WALTER  KRCLL  ») 
(*  CF  THE  ENERGISE  SCIENCE  DEPARTMENT  FUNDED  BY  THE  ARMY  MEDICAL  »> 
<*  RESEARCH  COMMAND  .  *) 


PROGRAM  FASTMOVt INPUT/, OUTPUT/, EMGFILE  ,  CUTFILE.TAPEO  ,  TAFE1  ) ? 


'  CONST 

TS  = 
TA8I1  =  ' 
TA8I2  -  ' 
TA8I3  =  ' 
TABI4  =  ' 
TA8I5  =  ' 
TA8IS  =  ' 
TA8I7  =  ' 
TA0I8  =  ' 
TABI9  =  ' 
TASI10=  ' 

TABI 11=  '  - 

TAB! 15=  ' 
TA8I 12=  ' 
TA8I13=  ' 

TAEI14=  '  TIME 
TA8I15=  '  (MSEC) 


FAST  MOVEMENT  PROJECT  ", 

DEPARTMENT  OF  EXERCISE  SCIENCE  '  ,' 

UNIVERSITY  CF  MASSACHUSETTS/AMHERST '  ? 

PRINCIPLE  INVESTIGATORS  I  DR.  WILLIAM  KILMER  AND  DR.  WALTER  KRCLL'? 
RESEARCH  ASSISTANTS  I  VALERIE  CCNGDCN  AND  JEAN  BOUCHER '  ? 

MODEL  I  E  /  G  /  PH.IDOT  <  EMG  /  TORQUE  /  ANGULAR  VELOCITY  )  '  ? 

MODEL  PARAMETERS  ?  '  ? 

TRAPEZOID  RULE  STEP  SIZE  I  '  ? 

RUNGE-KUTTA  RULE  STEP  SIZE!  '  ? 

A  K  3ETA1  BETA2  XI  MU  Q(AITCH)  '  ? 


DATE 

SUBJECT  ID  I 

EMG  EMG  PHI 
(BICEP) (TRICEP)  (POS) 


SUBJECT  AM/PM  SESSION  BLOCK  TRIAL 


PHI  PHI  '  PSIV  PS  I  IE  QE 
(VEL)  (ACC)  '  ; 


ps:  if 


TYPE 

WHICHCNE  =  (  FLEXCR  ,  EXTENSOR  )  ? 

REAL2D  =  ARRAY! 0..25,!.. 23  CF  REAL  ? 

SAVE IT  =  ARRAY! 1 . .70,0. .11  CF  REAL  ?  <*  NOTE!  ?J  =  MU  +  XI  +  DELAY (=20)  *! 


VAR  EHG3ICEP  I  PEAL2D  ", 

emgtricep:  realzd  : 
savephis  :  saveit  ; 


ET  ,  EB  .  PSIVE  ,  PSIIE  ,  PEI  IF  , 

QUEE  ,  QUEF  , 

PHI  ,  PHIDOT  ,  PHIDDQT  I  REAL  ? 

BETA1F  ,  BETA1E  , 

BETA2F  ,  EETA2E  , 

XIF  ,  XIE 
MUF  ,  MU&  • 

AF  ,  AE  ,  KF  ,  KE  . 

BETA  ,  CZETA1  ,  CZETA2  ,  CZETA  , 

GF  ,  QE  I  REAL  ? 

LASTPHI  ,  LASTPHIDCT  ,  LASTPHIDDOT  , 

PI  ,  P2  .  P3  .  CCF'AU  .  ”  .  MAN  .  DELAY  I  REAL  ? 

BACK  INTIME  .  I,MM,DD, YY,SU3,AP,SEa,cL,TR  I  INTEGER  ? 

HTRAP  ,  HRUNGE  ",  REAL  ",  U  H  FCR  TRAPEZOID  AND  RUNGE-KUTTA  RULES  •>) 
EMGFILE  ,  OUTFILE  ,  TAPEO  ,  TAFE1  I  TEXT  ? 

DAMPSET  I  BOOLEAN  ", 


ocedures  a n 


ion.; 


un ; 


procedure  reademcfile  : 

Reads  in  data  from  input  file  EU3FILE  and  daces  id 
in  Che  appropriate  array.  See  Input  ana  Output  secti 


FUNCTION  3ETEMG  : 

For  a  particular  time  this  function  calculates  the 
Smoothed  EMG  value  for  either  the  flexor  or  extensor 


FUNCTION  EMC  : 

Computes  Augmented  EMC  values. 


PROCEDURE  PUSH  : 

Saves  PHI , PH I DOT, and  PHIDDOT 


values  back  in  time. 


FUNCTION  CJE 


For  e r.uh  set 
the  'q1  used 


of  beta  values  this  function  calculates 
in  tha  hit-tau)  fur.cti on .  ( See  text). 


FUNCTION  AITCH  : 

The  reflex  and  mechanical  response 
extensor  is  returned  from  tnis  fun: 


f  the  flexor  or 
ion.  (.  HCt-teu)) 


FUNCTION  TRAPEZOID  : 

A  numerical  method  for  approximating  the  integral  wh 
determines  tne  value  of  PHI  from  FHIDOT  and  PHICOT 
from  PHIDDOT. 


FUNCTION  PS 17  : 

The  Force-Velocity  relation  is  returned  wnen  given 
PHIDOT  as  a  parameter. 


FUNCTION  P3II  : 

The  Integral  Factor  of  Torque  is  approximated  using 
Runge-ivutta  rnetnod  for  integration. 


iN .  .  ALL  I 

Initialisation  of  ell  constants  ,  variables  and  arra 
In  addition  the  headers  of  the  output  files  are  prin 
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PROCEDURE  READEMGFILE(VAR  EMG3ICZ?  :  REAL2D  I 

var  ekgtrioep:  realzd  ; 

VAR  NAXENGTINE  REAL  )  ' 

VAR  INDEX. TIKE. V ALU  INTEGER 

ENDFILE  I  3COLEAN 

BEGIN 

RESET (ENGFILE)  ,' 

index  :=  o  ; 
tine  ; 
valu  :=  2  ; 

ENDFILE  I-  FALSE 

READ (ENGFILE.  NM.DD. YY.SU3.AP.SES.BL. TR  ,  MAXENGTINE)  ! 

WHILE  NOT  ENDFILE  DO 
BEGIN 

READ (ENGFILE, 

EMGBICEPC INDEX , TIME] ,  EMGBICEPC INDEX, VALU] , 
EMGTRICEPC INDEX, TIME] , EMGTRICEPC INDEX, VALU] )  ! 

IF  EMGBICEPC INDEX. TIME]  >=  NAXENGTINE 
THEN  BEGIN 

ENDFILE  1=  TRUE  I 
INDEX  1=  INDEX  -  1 
END 

ELSE 

INDEX  1=  INDEX  +  1 
END  ; 

END,"  (*  END  OF  PROCEDURE  READENGFILE  »> 


FUNCTION  GETENG!  ATTIXE  ,  NAXENGTINE  I  REAL  !  VAR  ENGA  I  REAL2D  )I  REAL,' 

(*  THIS  FINCTIGN  RETURNS  THE  VALUE  OF  THE  ENG  AT  TINE  =  ATTINE 
USING  THE  FOLLOWING  FCRNULA  I 

F(X)  =  ((  F(B)  -  F(A)  )/(8-A))*(X-A)  +  F(A)  *> 

(*  THE  ENGA  HERE  IS  THE  ARRAY  HOLDING  THE  INTERVAL  VALUES  E.G.  ENG3ICEP  *) 

VAR  PREVICUSTINE  ,  “REVICUSVALU  I  REAL  ,' 

TIME  ,  VALU  ,  INDEX  I  INTEGER  ! 

BEGIN 

IF  ATTINE  <0.0 
THEN  GETENG  O.C 
ELSE  BEGIN 

index  ; =  i  ; 
tine  :=  i  ; 
valu  :=  2  ; 

(*  INCREMENT  INDEX  UNTIL  UPFER  BCUND  OF  INTERVAL  IS  FOUND  »> 

WHILE  ((  ATTINE  >  EMGAC INDEX, TINE!)  AND 
(  ATTINE  <=  NAXENGTINE)) 

DO  INDEX  :=  INDEX  +  1  ! 

(*  GET  LOWER  SOUND  OF  INTERVAL  *) 

IF  INDEX  )=  i 
THEN  BEGIN 

PREVIOUSTIKE  ENCAD  INDEX-1,  TINE]  ; 

PREVICUSVALUI =  EMGAC INDEX-1 .VALU] 

end; 

<♦  ATTINE  IS  EGUAL  TO  UPPER  SOUND  OR  WHOLE  INTERVAL  IS  THE  SANE  VALUE  *) 
IF  ((  ATTINE  *  EMGAC INDEX, TIME]  )  GR  (  EMGAC INDEX, VALU]  =  PREVICUSVALU  )! 
THEN 

GETENG  EMGAC  INDEX,  VALU] 

ELSE 

<«  UPPER  BCUND  IS  ZERO  THEREFORE  THE  WHOLE  INTERVAL  IS  ZEROES  *) 

IF  EMGAC INDEX, VALU]  =  0.0 
THEN  GETENG  ;=  O.C 
ELSE 

(»  ATTINE  IS  SOMEWHERE  INSIDE  THE  INTERVAL  *) 

GETENG  I =  <  <  EMGAC INDEX, VALUI-  PREVICUSVALU) 

/(  emgacindex.time;  -  previoustime) : 

*  (ATTINE  -  FREVIOUSTIME)  +  PREVICUSVALU 
END 

END; (5  END  OF  FUNCTION  GETENG  *) 
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FUNCTION  ENG!  T  .  TAU  I  REAL  ,*  NUSCLE  I  WHICH  ONE  )  I  REAL  I 

<»  RETURNS  THE  AUGMENTED  ENG  VALUE  FOR  EITHER  THE  FLEXOR  OR  EXTENSOR  *) 

VAR  BACK  I  INTEGER  ,' 

EMGA  I  REAL2S  1  ZERO  ,  ONE  *  TWO  .  X  I  REAL  ! 

BEGIN 

BACK  I’  TRUNC!  T  -  TAU  ♦  DELAY  )  ' 

ZERO  :=  SAVEPHISCSACX  >011 
ONE  I*  SAVEPHISI3ACX  .13.* 

TWO  .*  =  SAVEPHISI3ACX  >231 

IF  MUSCLE  =  FLEXOR 
THEN  BEGIN 

EMGA  1=  EMG8ICE?  1 

X  1=  ZERO  +  P2  *  ONE  ♦  P3  *  TWO 
END 

ELSE  BEGIN 

EMGA  1=  EMGTRICEP  1 

X  1=  -1.0  #  ZERO  -  P 2  *  ONE  -  P3  *  TWO 
END  ; 

IF  ((  ZERO  >=  0.0  )  AND  (GGFTAU  =  0.0)) 

THEN 

EMG  .'=  GETEMG  ( TAU »  MAX .  EMGA ) 

ELSE  BEGIN 

GCFTAU  :=  1.0  ; 

IF  X  >  0.0  (♦  CUSPED  BRACKETS  RULE  ») 

THEN 

EMG  .*=  (  PI  #  X  *  GOFTAU)  +  GETEMG ( TAU > MAX , EMGA ) 

ELSE 

EMG  Is  GETEMG (TAU, MAX, EMGA) 

END 

END  ;<♦  END  OF  FUNCTION  EMG  *) 


PROCEDURE  PUSH (  VAR  V0,V1,V2  I  REAL  1  VAR  ONTO  I  SAVEIT)  1 
(*  PLACES  THE  VALUES  V0,V1,V2  IN  A  STACK  NAMED 

ONTO  WHICH  SAVES  UPTO  70(MU+XI+20)  VALUES  BACK  IN  TIME  *) 

VAR  I  I  INTEGER  1 

BEGIN 

FCR  I  :=  0  TO  BACK INTIME  -  2 
DO  BEGIN 

CNT0C3ACXINT:ME-I,03  :=  CNT0C3ACXINTIME-I-1,03  ; 
ONTOIBACKINTIME-1 , 1 3  CNTOCEACX INTIME-I-1 , 13  1 
ONTO  C  SACK  INTI  ME- 1 ,  Z  3  1=  CNT0C3ACX INTIME-I-1 ,23  ’,  ' 
END  ; 

ONTO! I ,03  :=  VO  ; 

onto:i,i:  :=  vi  ; 

ONTO Cl, 2 3  I*  V2 

END  ;<*  END  OF  PROCEDURE  PUSH  *) 


FUNCTION  GUEt  MUSCLE  I  HHICHGNE  )  l  SEAL  J 

(*  THIS  FUNCTION  FINES  THE  VALUE  OF  CUE  SUCH  THAT  THE  AITCH  FUNCTION 
CONDITIONS  ARE  SATISFIED  IE.  GUE  TIMES  THE  INTEGRAL  =1.0  * 

VAR 

sigma  ,  bi  ,  bo  ,  ::i  :  real  ; 


BEGIN 

IF  MUSCLE  =  FLEXOR 
THEN  BEGIN 

Bl  I*  SETA1F  J 
B2  I*  SETA2F  < 

XI  1=  XIF  ; 

END 

ELSE  BEGIN 

Bl  .'=  BETA1E  ; 

B2  :=  BETA2E  » 

xi  :=  xie  ; 

END  ; 

sigma  :=  xi  ; 

GUE  :=  <  1.0  /  <((  EXP(  -1.0*B2*SIGMA)  -  1.0)  /  e2  )  ♦ 

<<  1.0  -  EXP(  -1.0*B1*SIGMA  )  >  /  Bl  ))  ) 

END  ; 

(*  END  OF  FUNCTION  GUE  *) 


FUNCTION  AITCH (  TIME  I  REAL  1  MUSCLE  I  WHICHONE  )I  REAL  I 
(*  Th'E  AITCH  FUNCTION  RETURNS  THE  REFLEX  AND  MECHANICAL 
RESPONSE  VALUE  FOR  TIME  =  T  -  TAU  •) 

VAR 

HH  .  G  ,  ei  ,  B2  .  XI  ,  MU  ,  TSTAR  I  REAL 
BEGIN 

IF  MUSCLE  =  FLEXOR 
THEN  BEGIN 

Q  .*=  OF  ! 

31  :=  BETA1F  ; 

B2  1=  3ETA2F  ", 

xi  :=  xif  ; 
mu  :=  muf  ; 

END 

ELSE  BEGIN 

Q  1=  GE  ; 

Bl  :=  BETA1E  ; 

82  1=  BETA2E  1 

xi  :=  xie  ; 

MU  :=  MUE 
END  ; 

IF  <(  TIME  <  MU  )  OR  (  TIME  >  (XI  ♦  MU)  )) 

THEN  AITCH  1=  0.0 
ELSE 
BEGIN 

TSTAR  :=  TIME  -  MU  ' 

HH  :=  G  *  (  EXP(  -1.0*81*TSTAR  )  -  EXP(-1 .0*B2*T5TAR  )  )  ! 
AITCH  :=  HH  ; 

(*writeln( '  H( ',time:4:o, ')  =  ',hh:io:5)  ;*) 

END 

end; 

(#  END  FUNCTION  AITCH  #) 


FUNCTION  TRAFEZQIDt  HTRAP  .  LASTCNE  .  LASTFUN  ,  FUN  I  REAL)  :  REAL  ’ 

(*  WE  ARE  USING  THE  TRAPEZOID  RULE  TO  GET  PHI  FROM  PHIDCT  FRCM  PHIDDCT 
E.G.  FH I (M)  =  PHKN-l)  +  .5»HTRAP*<  PHIDCT:N-1)  +  PHIDQT(N)  ) 
LASTCNE  LASTFUN  FUN 


BEGIN 

TR.APEZGID  .*=  LASTCNE  +  (  0.5*  HTRAP  *  (  LASTFUN  +  FUN  )  ) 
END  ! 

(*  END  OF  FUNCTION  TRAPEZOID  *) 


»  ... 


FUNCTION  PSIVt  DOT  I  REAL  )  I  REAL  ", 

(*  P SIV  .  THE  FORCE  -  VELOCITY  RELATIONSHIP 
(»  PSIV  =  PSIVE  =  PSIVF 
VAR  X  :  REAL  ", 

BEGIN 

X  :=  1.0  -  (  DOT  /  1000.0  )  ", 

IF  X  >=  0.0 
THEN  PSIV  1=  X 
ELSE  PSIV  1=  0.0 

end; 

<#  END  OF  FUNCTION  PSIV  *> 


FUNCTION  PSIIC  MUSCLE  I  WHICHONE  ",  T  I  REAL)  I  REAL 
(*  THE  RUNGE-XUTTA  METHOD  FOR  INTEGRATION  IS  USED  HERE 
REF I  NUMERICAL  METHODS  ;  DAHLGUIST  ET.AL.  .PRENTICE 
HALL,  1974.  PP  34S-7  *) 


VAR  INTEGRAL  ,  X  ,  KCNE  ,  KTWO  ,  K THREE  >  A  •  H  , 
XI  ,  MU  I  REAL 
I  ,  NUMINTS  I  INTEGER  ", 


FUNCTION  F(  TAU  ,  T  I  REAL  ",  MUSCLE  I  WHICHONE  )  I  REAL 
VAR  EMGA  I  REALOD  ,*  K  I  REAL  ", 

BEGIN 

IF  MUSCLE  =  FLEXOR 
THEN  K  1=  KF 
ELSE  K  1=  KE  ; 

F  1=  K  *  EMG(  T  ,  TAU  ,  MUSCLE  ) 

*  AITCH(  T-TAU  .  MUSCLE  ) 

END  J  t#  END  OF  FUNCTION  F  IN  PS! I  *) 


BEGIN  (*  PSII  MAIN  SECTION  *) 

IF  MUSCLE  =  FLEXOR 
THEN  BEGIN 

XI  1=  XIF  ; 

mu  :=>  muf  ; 

A  :=  AF 
END 

ELSE  BEGIN 

xi  :=  xie  ; 
mu  :=  mue  ; 

A  AE 
END  ; 

H  :=  HRUNGE  ; 

integral  :=  o.o  ; 

NUMINTS  :=  TRUNCC  (XI/H)  -  1.0  )  ! 

FOR  I  I-  0  TO  NUMINTS 
DO  3EGIN 

X  :»  T  -  MU  -  XI  +  (H  *  I)  ;  (*  TAU  IN  TEXT  *) 

KONE  :=  Ft  X  .  T  .  MUSCLE  )  ", 

KTWO  .*=  Ft  X  +  H/2.0  ,  T  ,  MUSCLE  ) 

KTHREE  1=  F(  X  +  H  ,  T  ,  MUSCLE  )  ! 

INTEGRAL  I*  INTEGRAL  +  (  (H/S.O)  *  (  KCNE  +■  4.0  ♦  KTWO  +  X THREE 
END  ", 

?SII  Is  A  »  CCS!  PH  I  *  0.C17iS,C''!  *  INTr"RAL 

“phi  in  radians"  2*p:  7“cgo  *) 

END  ;  (♦  END  OF  FUNCTION  PSII  ♦) 


PROCEDURE  INITALL  J 


44 


I  INTEGER  ; 

BEGIN 

DAHPSET  FALSE  I 

reurite(Outfile)  ; 

REWRITE  IT APEC)  i 
REwRITE( TAPED  ! 

CZETA1  :=  0.01  ; 

CZETA2  :=  0.09  1 

BETA1E  1=  0.10  "  8ETA2 £  .*=  0.15  ,' 

SETA1F  :=  0.10  ;  BETA2F  0.15  1 

xie  :=  4o.o  ; 

XI F  :=  40.0  ,' 

mue  :=  go.o  : 
muf  :=  so.o  ; 

GE  :=  QUE<  EXTENSOR  )  ; 

OF  :*  GUE(  FLEXOR  )  ! 

htrap  :=  i.o  ; 

HRUNGE  :=  0.5  ; 

pi  :*  o.oo  ; 

P2  :=  o.o  ; 

P3  :=  oo.o  ; 

GCFTAU  :=  0.0  ; 

ke  :=  0.40  ; 
kf  o.65  ; 
ae  :=  o.io  ; 
af  :=  o.i  ; 


phi  :=  so.o  ; 
ph i dot  :=  o.o  ; 
phiddct  :=  o.o  ; 
lastphi  :=  so.o  ; 

LASTPHIDGT  1=  0.0  1 
LASTPHIDDOT  1=  0.0  ! 

DELAY  D  20.0  ! 

3ACXINTIXE  TRUNC (  HUE  +  XIE  +  DELAY  1  I 

FOR  I  1=  1  TO  BACK  INTINE 
DO  BEGIN 

SAOEPHISCI  .01  :=  0.0  <*  PHI  *)  ! 

saoephisidd  :=  o.o  u  phidot  *>  ; 

SAVEPHISCI.2I  1=  0.0  <*  PHIDDGT  *) 
END  ; 


READENCFILE!  ENGBI 


ENGTRICE?  ,  NAX  I 


NR  I TELN  <  3UTF  I LE )  ’  UR I TELN (CUTFILE)  ! 
uritzln<cutfile,  tabid  ;  ur:teln(CUTf:lz)  ; 

UR I TELN (CUTFILE  ,  TA3X2) 

URITELN(CUTFILE  ,  TABI31  1  UR I TELN (CUTFILE! 

URITELN(CUTFILE  ,  TAS I 4 )  !  UR I TELN ( CUTF I LE  -  TA3I5)  ! 

WRITELN(CUTFILE)  URITELN ( CUTFILE ) 

UR I TELN (CUTF I LE  .  TA3I3)  )  UR I TELN ( CUTr I LE )  !  URITELN  ( CUTFILE  1  ! 

UR 1 TELN ( CUTF ILE  ,  TA3I7, '  DELAY  1  '.DELAY! 

wRI TELN  1 CUTFILE , TAB  IS , HTRAP ! 4 ! 1 . '  PI  ! ' rPi:S!3» '  P 2  !  ',? 2!3.u, 

'  P3  :  ',?3:g:3  >  ; 

X3 1  TELN  <  Du  TF I  LE  .  TA3IS  ,  HRUNGE.M.'D'  CZETA1  !  '  .DZETA1  !3.‘3. '  CZETA2 
WR I TELN ( OUTF I LE )  ! 

UR I TELN (CUTFILE  -TABI10)  1 

uritcL.v(gutfile.  'EXTENSOR  '  ,ae:4:2,xe:s:2,3etaie:5:z, 

3ETA2EI5I2,  XIE!3!2,  NUEISIZ.  3E:S!5  )  I 
UR I TELN ( CUTFILE  # 'cLEXCR  '  ,AF ! 4 ! Z » KF ! 3 !2 , 3ETA1 F ! 5! C . 

3ETA2F:5!2,  X:F!3!2,  NUFISIZ,  3F:8!5  )  ! 

UR I TELN (CUTFILE)  ! 

UR  I  TELN  (CUTFILE,  TA3IID  !  UR  I  TELN  (CUTFILE)  ! 

UR  I  TELN  ( DUTF ILE )  I 

ur;teln(Cutfile,tasii5)  ; 

URI  TELN  (CUTFILE,  TAoIlZ  ,  ,r.r:2,TS,DD:2,T3,  YY!2,  '  SUSIE,  ' 

'  ',3L!2,  '  ' , TR !3)  I  UR I TELN (CUTFILE)  ’ 

ur;teln(dl"file  ,  tasii3)  : 

UR  I  TELN  ( CUTF I  LE  ,  TA3H4)  I 
UR  I  "ELN!  CUTFILE  ,  TA3I15)  '• 
wRITiLN'CUTF iLi.)  ,  URITELNvCUirILc, »A3I*1)  J 
END  ,'  (*  END  CF  PROCEDURE  INI  TALL  +) 


',CZETA2!5 


',a?:d  ' 


URITELV(CUTFILE)  ! 


•C»l 


THE  MAIN  PROGRAM  _ 


(*. 


BEGIN 

INITALL  ! 

t  o.o  ; 

WHILE  T  <=  MAX 
00  BEGIN 

ET  1=  GETEMGfT,  MAX,  EMGTRICEP)  ! 

EB  1=  GETEMG  ( T ,  MAX,  EMGB I  CEP) 

psive  :=  fsiv<  PHioar  )  ; 

PSIIE  :=  PSII!  EXTENSOR  ,T  )  ! 

P SI IF  .'  =  PSII (  FLEXOR  ,T  )  J 
QL'EE  :=  PSIVE  *  PSIIE  J 

quef  :=  psive  *  psiif  ; 

IF  ((  PHI  >  0.0  )  AND  (  PHIDDT  >0.0  )) 

THEN  DAMPSET  :=  TRUE  ! 

IF  DAMPSET 

THEN  C2ETA  .*=  C2ETA2 
ELSE  CZETA  :=  CZETA1  ‘ 

PHIDDOT  1=  GUEE  -  GUEF  -  CZETA*PHIDOT  ; 

(*  ZETA  *) 

PHIDOT  1=  TRAPEZOID!  HTRAP  ,  LASTPHIDOT  ,  LASTPHIDDOT 

PHI  1=  TRAPEZOID!  HTRAP  ,  LASTPHI  ,  LASTPHIDOT 

!*  WRITE  TO  OUTPUT  FILE*) 

IF  !  TRUNC(T)  MOD  10  )  =  0 
THEN  BEGIN 

writeln(Tapeo,E3:b:4,et:b:4  j  ; 
writeln(tapei,phi:3:4,phibot:b:4  )  ; 
writeln(outfile,t:s:o,  '  ', £8:7:3, £T:a.'3,PHr:s:4, 

PHiD0T:3:4,PHi3DGT:s:4,psiv£:3:4,PsiiE:s:4,auEE:s:4 

psiif:9:4,guef:s:4)  end  ; 

PUSH!  PHI, PHIDOT, PHIDDOT, SAVEPHIS)  ! 

IF  T  =  555 

THEN  FOR  I  1=  i  TO  BACK  INTIME 
DO 

WRITELN! 1 14, SAOEPHISCI, 0i:8:3, SAVEPHISCI, 12:9:3, 

savephisci, 22:9:3)  ! 

LASTPHI  :=*  PHI  ; 

LASTPHIDOT  :=  PHIDOT  ! 

LASTPHIDDOT  :=  PHIDDOT 

t  :=  t  +  t.o  ; 

END  ; 

END. 


PHIDDOT  ) 
PHIDOT  ) 


EMGFILE  :  Input  file  contains  the  Smoothed  EMG  data 


OUTFILE  :  Output  file  showing  values  used/computed  in  table  form 
TIME  ,  ET  ,  EB  ,  PHI  ,  PHIDOT  ,  PHIDDOT  ,  PSIV 
PSIIE  ,  QE  ,  PS I IF  ,  OF  . 


TAPEO  :  Output  file  holding  values  of  ET  and  E3  for  PLOT  Program 


TAPE1  :  Output  file  holding  values  of  PHI  and  PHIDOT  for 
PLOTPHI  Program. 


PLOT  Program  :  Written  in  FORTRAN  this  Program  plots  the 

EMG  values  over  time.  TAPEO  is  its  input  file. 


PLOTPHI  Program  :  Written  in  FORTRAN  this  Program  plots  the 
Position  and  velocity  values  over  time. 
TAPE1  i3  its  input  file. 


EMGFILE 


Line  number  1  contains  the  Subject  Identificaion  numbers 
and  the  maximum  time  of  trial. 

Linas  2  through  (a  possible)26  hold  the  following  data  : 

Column  1  :  The  Biceps  stop  time  of  interval 

Column  2  :  The  Biceps  value  at  end  of  interval 

Column  3  :  The  Triceps  stop  time  of  interval 

Column  4  :  The  Triceps  value  at  end  of  interval 

Note  :  If  either  the  Biceps  or  Triceps  has  fewer  intervals  the  table 

Must  be  padded  so  that  each  column  has  the  same  number  of  entries. 

The  SUBJECT  ID.  is  made  up  of  eight(8)  variables  : 

Mounth, Day, Year, Subject  number, AM  or  PM, Session, 31ock, and  Trial  numbers. 

EXAMPLE  : 

On  November  14,1980  subject  41  exercised  in  the  AM  for 
the  fourth  day  first  session,  bloi'x  one  and  15th  trial  with 
the  following  results  after  500  milliseconds  : 


14  80 

10  4  1 

15  500. 

0 

0.0 

0.0 

0.0 

0.0 

64.0 

0.65 

166.0 

0.07 

96.0 

0.65 

136.0 

0.58 

323.0 

0.0 

228.0 

0.95 

324.0 

0.09 

3  00.0 

0.95 

350.0 

0.13 

350.0 

0.0 

500.0 

0.0 

500.0 

0.0 

OUTFILE 


Contains  the  parameters  used  and  the  values  computed  in  PROGRAM  FASTMCV. 


Sample  Header 


FAST  MOVEMENT  PROJECT 

DEPARTMENT  OF  EXERCISE  SCIENCE 
UNIVERSITY  OF  MASSACHUSETTS/AMHERST 

PRINCIPLE  INVESTIGATORS  I  DR.  WILLIAM  KILMER  AND  DR.  WALTER  KRCLL 
RESEARCH  ASSISTANTS  I  VALERIE  CCNGDON  AND  JEAN  BOUCHER 


MODEL  I  E  /  Q  /  PHIDOT  (  EMG  /  TORQUE  /  ANGULAR  VELOCITY  ) 


MODEL  PARAMETERS  1  DELAY  I  20.00 

TRAPEZOID  RULE  STEP  SIZE  I  1.0  PI  .'  0.030  P2  I  S.COO  P3  I  2 0.000 

RUNGE-KUTTA  RULE  STEP  SIZE.*  0.5  CZETA1  I  0.010  CZETA2  1  0.500 

A  K  BETA1  BETA2  XI  MU  Q(AITCH) 

EXTENSOR  0.10  0.40  0.09  0.15  40.00  SO. 00  0.24054 
FLEXOR  0.10  0.65  0.09  0.15  40.00  60.00  0.24054 


DATE  SUBJECT  AM/PM  SESSION  BLOCK  TRIAL 
SUBJECT  ID  1  11/14/80  104  1  15 


TIME  EMG  EMG  PHI  PHI  PHI  PSIV  PSIIE  QE  PSIIF 

(MSEC)  (BICEP) (TRICEP)  (PCS)  (VEL)  (ACC) 


TAPEO  : 


Column  1  holds  the  value  of  the  Smoothed  EMG  (BICEPS).  (XXX.XXXX) 
Column  2  holds  the  value  of  the  Smoothed  EMG  (TRICEPS).  (XXX.XXXX) 


TAPE!  I 


Column  1  holds  the  value  of  PHI(t).  (XXX.XXXX) 

Column  2  holds  the  value  of  PHIDOT(t).  (XXX.XXXX) 


PLOT 
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PROGRAM  PLOT< INPUT. OUTPUT, TAPEO)  _  „  „ 

#  PLOTS  TWO  GRAPHS  SIMULTANEOUSLY  WITHIN  THE  RANGE  -.5  TO  1.5 

*  TOTAL  OF  51  POINTS. 

DIMENSION  DATA<51,3.3> 

PRINT  20 

20  FORMAT  (//>"  MODEL  I  E  /  Q  /  PHIDOT  (  JMG  /  TORQUE/// 

+  "ANGULAR  VELOCITY  )*,//,“  E  .  EXTENSOR  < TRICEPS )  ./» 

+■  F  I  FLEXOR  (BICEPS)  “»/."  1  1  INDICATES  THE  1.0  LINE  » 

+  /,"  0  I  INDICATES  THE  0.0  LINE  ‘.//) 

DO  10  Isl>51 
READ! 0.50)  DAT1  ,  DAT2 
IF(E0F(0))99.57 

57  IF(DAT1.GT.1.5)  DAT1  =  1.5 

IF(DAT1.LT.~,5)  DAT1  *  -.5 
56  IFIDAT2.GT.1.5)  DAT2  *  1.5 

IFtDAT2.LT.-. 5)  DAT2  =  -.5 
I ALT  *  10 

DATA(I.l.l)  *  FLOAT! I*IALT)  -  10.0 

DATA! I, 1.2)  =  DAT1 

DATA ( 1.1.3)  =  IHF 

DATA< 1.2.2)  *  DAT2 

DATA(I.2.3)  =  1HE 

DATA( I .3,2)  *  1. 

DATA( 1.3.3)  =  1H1 
10  CONTINUE 

CALL  PLQTIN(2AT- .-51 ,51, 3, -.5, 1.5, 1H0.75.0) 

50  FORMAT (2FB. 4) 

go  CMn 

SUBROUT INE  PLOT  I N  ( D  AT  A ,  NX ,  M  AXX ,  N  YSETS .  YM I N ,  YMAX ,  AX  I S ,  I PW ,  I RCL  > 
DIMENSION  DATA(MAXX.NYSETS.3),PL0T2(120) 

DO  B  L=1 , 120 
3  PLQT2(L)=1H 

IF(YMAX.G£.YMIN)30T010 

2  FORMATO^'mIn' y'^IS  GREATER  THAN  MAX  Y*/*  MIN  Y*  *,G17.10,*  MAX  Y= 
X  G17.10) 

RETURN 

10  Y1NC* tYMAX-YMIN) /FLOAT! IPW-1Z) 

DO  100  JX=1 .NX 
DO  20  JY=1 .NYSETS 
Y=DATA( JX, JY.2) 

IF ( Y . GE . YM I N . AND . Y . LE . YMAX ) G0T030 

3  FCRMATu\  VALUE  DOES  NOT  FIT  WITHIN  DEFINED  Y  RANGE*. G17. 10) 

RETURN 

30  INDX=INT( .5+<Y-YMIN)/YINC) 

IF( INDX.LT . 1 ) INDX*1 
IF(INDX.GT.120)INDX=120 
PL0T2 ( INDX ) =DATA ( JX , J Y , 3 ) 

20  CONTINUE 
C 

IF(AXIS.Efl.0)G0T040 

IF(YMIN.GT.0)G0T040 

IXAXIS*(0-YMIN)/YINC 

IF(PL0T2(IXAXIS) .E3.1H  )PL0T2(IXAXIS)*AXIS 
40  L£N=LENGTH ( PLQT2 ) 

PRINT  5,DATA( JX, 1.1), (PL0T2(K) ,K=1 ,LEN) 

5  FORMAT (Gil .5, IX. 120A1 ) 

DO  50  L=  1.120 
50  PLQT2(L)*1H 
100  CONTINUE 
RETURN 
END 

INTEGER  FUNCTION  LENGTH! I ARRAY) 

C  RETURNS  LENGTH  CF  NON-BLANK  PART  OF  IARRAY. 

DIMENSION  IARRAY (120) 

LENGTH *0 
DO  100  1*1,120 

IF! IARRAY! J). NS. 1H  )  GOTO  200 
100  CONTINUE 
200  LENGTH* J 
RETURN 
END 
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PROGRAM  PLOTPHI ( INPUT , OUTPUT/ TAPE1 ) 

*  PLOTS  TWO  GRAPhS  SIMULTANEOUSLY  WITHIN  THE  RANGE  -150.  TO  150. 

♦  TOTAL  OF  50  POINTS. 

DIMENSION  DATA151 ,3,3) 

PRINT  20 

20  FORMAT!//, “  MODEL  I  E  /  0  /  PH I DOT  (  EMG  /  TORQUE  / 

+  “ANGULAR  VELOCITY  )",//,"  PI  POSITION  (PHI)  “,/, 

+  V  I  VELOCITY  *  100  (PHIDOT  *  100)  ", 

+  /,"  .  I  INDICATES  SO  DEGREE  LINE  ",/, 

‘  00000  I  INDICATES  0  DEGREE  LINE  “,//) 


SO 


DO  10  1*1,51 
READ< 1 ,50)  DAT1  ,  DAT2 
IF(ECF(1))99,57 
57  DAT2  *  DAT2  *  100.0 

IF1DAT1 .GT. 150. )  DAT1  =  150. 
IF(DAT1.LT.-150.)  DAT1  *  -150. 

56  IF(DAT2.GT.1S0. )  DAT2  *  150. 

IF(DAT2.LT.-150. )  DAT2  *  -150. 

I ALT  *  10 

FLOAT! I* I ALT)  -  10.0 
DAT1 
1HP 
DAT2 
1HV 
SO. 

1HI 


DATA! I, 1,1) 

DATA! I, 1,2) 

DATA! I, 1,3) 

DATA! I ,2,2) 

DATA! I ,2,3) 

DATA! I, 3, 2) 

DATA! 1,3,3) 

10  CONTINUE 

CALL  PLOT IN (DATA, 51 ,51 ,3,-150. , 150. , 1H0,75,0) 

50  FORMAT (2F9. 4) 

99  END 

SUBROUTINE  PLQTIN! DATA . NX , HAXX , NYSETS , YMIN , YMAX , AXIS , I PW , IRCL ) 
DIMENS ION  D ATA ( MANX , NYSETS , 3 ) , PL0T2 (120) 

DO  9  L* 1,120 
9  PL0T2!L)=1H 

IF ( YMAX. GE. YMIN) GCT010 
PRINT  2, YMIN, YMAX 

2  FORMAT ( *  MIN  Y  IS  GREATER  THAN  MAX  Y*/*  MIN  Y=  *,G17.10,*  MAX  Y* 
X  G17.10) 

RETURN 

10  Y I NC  = ( YMAX- YM I N ) /FLOAT ( 1 PW- 1 2 ) 

DO  100  JX=1 ,NX 
DO  20  JY=1, NYSETS 
Y=DATA(JX, JY,2) 

IF(Y. GE. YMIN. AND. Y.LE. YMAX)G0T030 
PRINT  3, Y 

3  FORMAT ( *  Y  VALUE  DOES  NOT  FIT  WITHIN  DEFINED  Y  RANGE*, Gl7. 10) 

RETURN 

30  INDX*INT(.5+(Y-YMIN)/YINC) 

IF(  I.NDX.LT .  1 )  INDX=1 
IF( INDX.GT. 120) INDX=120 
PL0T2 ( INDX) *DATA ( JX , J Y , 3) 

20  CONTINUE 

IF ( AX IS . EG . 0 ) G0T040 
IF(YMIN.GT.0)GQT04O 
IXAXIS*(0-YMIN) /YINC 

IF(PL0T2( IXAXIS) .EG. IH  )PL0T2( IXAXIS)*AXIS 
40  LEN=L£NGTH( PLQT2) 

PRINT  5,2ATA( JX, 1 , 1 ) , (PLQT2(K) ,K=1,L£N) 

5  F0RMAT(G11 .5, IX, 120A1 ) 

DO  50  1-1,120 
50  PL0T2(L)*1H 

100  CONTINUE 
RETURN 
END 

INTEGER  FUNCTION  LENGTH ( I ARRAY) 

RETURNS  LENGTH  CF  NCN-3LANK  PART  CF  I ARRAY. 

DIMENSION  IARRAYI 120) 

L£NGTH*0 
DO  100  1*1,120 
j- I 

IF(IARRAY(J).NE.1H  )  GOTO  200 
100  CONTINUE 
2 00  LENGTH* J 
RETURN 
END 
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